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Abstract 

Using tools from algebraic geometry and Grobner basis theory we solve 
two problems in network coding. First we present a method to determine 
the smallest field size for which linear network coding is feasible. Second 
we derive improved estimates on the success probability of random lin- 
ear network coding. These estimates take into account which monomials 
occur in the support of the determinant of the product of Edmonds ma- 
trices. Therefore we finally investigate which monomials can occur in the 
determinant of the Edmonds matrix. 

Keywords. Distributed networking, linear network coding, multicast, 
network coding, random network coding. 

1 Introduction 

In a traditional data network, an intermediate node only forwards data and 
never modifies them. Ahlswede et al. [T] showed that if we allow intermedi- 
ate nodes to process their incoming data and output modified versions of them 
then maximum throughput can increase, and they also showed that the maxi- 
mum throughput is given by the minimum of maxflows between the source node 
and a sink node for single source multicast on an acyclic directional network. 
Such processing is called network coding. Li et al. [TH] showed that computa- 
tion of linear combinations over a finite field by intermediate nodes is enough 
for achieving the maximum throughput. Network coding only involving linear 
combinations is called linear network coding. The acyclic assumption was later 
removed by Koetter and Medard [5]. 

In this paper we shall concentrate on the error-free, delay-free multisource 
multicast network connection problem where the sources are uncorrelated. How- 
ever, the proposed methods described can be generalized to deal with delays as 
in [7] . The only exception is the description in Section [7] 

Considering multicast, it is important to decide whether or not all receivers 
(called sinks) can recover all the transmitted information from the senders 
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(called sources). It is also important to decide the minimum size q of the finite 
field F 9 required for linear network coding. 

Before using linear network coding we have to decide coefficients in linear 
combinations computed by intermediate nodes. When the size q of a finite field 
is large, it is shown that random choice of coefficients allows all sinks to recover 
the original transmitted information with high probability [7j. Such a method 
is called random linear network coding and the probability is called success 
probability. As to random linear network coding the estimation or determination 
of the success probability is very important. Ho et al. [7J gave a lower bound 
on the success probability. 

In their paper [9], Koetter and Medard introduced an algebraic geometric 
point view on network coding. As explained in [3 , computational problems 
in algebraic geometry can often be solved by Grobner bases. In this paper, 
we shall show that the exact computation of the minimum q can be made by 
applying the division algorithm for multivariate polynomials, and we will show 
that improved estimates for the success probability can be found by applying 
the footprint bound from Grobner basis theory. These results introduce a new 
approach to network coding study. As the improved estimates take into account 
which monomials occur in the support of the determinant of a certain matrix 
[7J we study this matrix in details at the end of the paper. 

2 Preliminary 

We can determine whether or not all sinks can recover all the transmitted infor- 
mation by the determinant of some matrix [7J. We shall review the definition 
of such determinant. Let G = (V,E) be an directed acyclic graph with possible 
parallel edges that represents the network topology. The set of source and sink 
nodes is denoted by S and T respectively. Assume that the source nodes S 
together get h symbols in F q per unit time and try to send them. 

Identify the edges in E with the integers 1, . . . , \E\. For an edge j — (it, v) 
we write head(j) = v and tail(j) = it. We define the \E\ x \E\ matrix F = (fij) 
where fo j is a variable if head(i) = tail(j) and fcj = otherwise. The variable 
fi j is the coding coefficient from % to j. 

Index h symbols in F q sent by S by 1, . . . , h. We also define an h x \E\ 
matrix A = (dij) where ajj is a variable if the edge j is an outgoing edge from 
the source s € S sending the z-th symbol and djj = otherwise. Variables Ojj 
represent how the source nodes send information to their outgoing edges. 

Let X(l) denote the Z-th symbol generated by the sources S, and let Y(j) de- 
note the information sent along edge j. The model is described by the following 
relation 

h 

2=1 i:head(i)— tail(j) 

For each sink t G T define an h x \E\ matrix B t whose entry bt.i,j is 
a variable if head(j) = t and equals otherwise. The index i refers to the i-th 
symbol sent by one of the sources. Thereby variables bt,ij represent how the 
sink t process the received data from its incoming edges. 

The sink t records the vector 
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where 

j:head(j)=t 

We now recall from [7] under which conditions all informations sent by the 
sources can always be recovered at all sinks. As in [7] we define the Edmonds 
matrix M t for t G T by 

m <={i-f bj)- W 

Define the polynomial P by 

p=ni M *i- ( 2 ) 

P is a multivariate polynomial in variables fij, ciij and &t,i,j- Assigning a value 
in F g to each variable corresponds to choosing a coding scheme. Plugging the 
assigned values into P gives an element k £ F g . The following theorem from [7] 
tells us when the coding scheme can be used to always recover the information 
generated at the sources S at all sinks in T. 

Theorem 1. Let the notation and the network coding model be as above. As- 
sume a coding scheme has been chosen by assigning values to the variables fij, 
aij and bt^i.j- Let k be the value found by plugging the assigned values into 
P. Every sink t £ T can recover from &W the informations X(l), . . . , X(h) no 
matter what they are, if and only if k ^ holds. 

Proof. See [7J. □ 



3 Computation of the Minimum Field Size 

We shall study computation of the minimum symbol size q. For this purpose we 
will need the division algorithm for multivariate polynomials [31 Sec. 2.3] to pro- 
duce the remainder of a polynomial F(X\ , . . . , X n ) modulo (Xf —X±,..., X% — 
X n ) (this remainder is independent of the choice of monomial ordering). We 
adapt the standard notation for the above remainder which is 

F(X U ...,X n ) rem (X« - X x , . . . , XI ~ X n ). 

The reader unfamiliar with the division algorithm can think of the above remain- 
der of F(X\, . . . , X n ) as the polynomial produced by the following procedure. 
As long as we can find an Xi such that Xf divides some term in the polynomial 
under consideration we replace the factor Xf with Xi wherever it occurs. The 
process continues until the X^-degree is less than q for all i = 1, . . . , n. It is 
clear that the above procedure can be efficiently implemented. 

Proposition 2. Let F(X\, . . . , X n ) be an n-variate polynomial overF q . There 
exists an n-tuple (x%, . . . , x n ) £ F™ such that F(xi, . . . , x n ) ^ if and only if 

F{X 1 , ...,X n ) rem ( Xf - X 1 , . . . , X« - X n ) ? 0. 

Proof. As a q = a for all a £ F q it holds that F(Xi, . . . ,X n ) evaluates to the 
same as . . . ,X n ) := F(X 1 ,...,X n ) rem (Xf - X u . . . , X*-X n ) in every 

(xi, . . . , x n ) € F™. If R(Xi, . . . , X n ) — therefore F(X±, . . . , X n ) evaluates 
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to zero for every choice of {x\, . . . , x n ) € F™. If R{X\, . . . , X n ) is nonzero we 
consider it first as a polynomial in F q {X\, . . . , X n -i)[X n ] (that is, a polynomial 
in one variable over the quotient field F q (Xi, . . . , X n _i)). But the X„-degree 
is at most q — 1 and therefore it has at most q — 1 zeros. We conclude that 
there exists an x n € F g such that i?(Ai, . . . , X n -i,x n ) £ F g [Ai, . . . , X„_i] is 
nonzero. Continuing this way we find (#1, . . . , x n ) such that R{x,\, . . . , x n ) and 
therefore also F{x-y, . . . , x„) is nonzero. □ 

From [7, Th. 2] we know that for all prime powers q greater than |T| linear 
network coding is possible. It is now straightforward to describe an algorithm 
that finds the smallest field F q of prescribed characteristic p for which linear 
network coding is feasible. We first reduce the polynomial P from ([2|) modulo 
the prime p. We observe that although P is a polynomial in all the variables 
b t ,i,j, fi.j the variable b ti ij appears at most in powers of 1. This is so 
as it appears at most in a single entry in M t and does not appear elsewhere. 
Therefore F q can be used for network coding if P rem p does not reduce to zero 
modulo the polynomials ajj — dij, ffj — fij. To decide the smallest field F q 
of characteristic p for which network coding is feasible we try first F g = F p . If 
this does not work we then try F p 2 and so on. To find an F q that works we 
need at most to try |_log_(|T|)J different fields as we know that linear network 
coding is possible whenever q > \T\. 

Note that once a field F q is found such that the network connection problem 
is feasible the last part of the proof of Proposition 1 describes a simple way of 
deciding coefficients (xi, ■ ■ . , x n ) G F™ that can be used for network coding. 

From [U Sec. 7.1.3] we know that it is an NP-hard problem to find the 
minimum field size for linear network coding. Our findings imply that it is 
NP-hard to find the polynomial P in 

4 Computation of the Success Probability of Ran- 
dom Linear Network Coding 

In random linear network coding we from the beginning fix for a collection 

KC{l,...,h}x{l,...,\E\} 
the djj's with € K and also we fix for a collection 

JC{l,...,|£|}x{l,...,|£|} 

the /ij's with (i, j) e J. This is done in a way such that there exists a solution to 
the network connection problem with the same values for these fixed coefficients. 
A priori of course we let a^j = if the edge j is not emerging from the source 
sending information j, and also a priori we of course let fi j = if j is not an 
adjacent downstream edge of i. Besides these a priori fixed values there may 
be good reasons for also fixing other coefficients dij and fij [7|. If for example 
there is only one upstream edge i adjacent to j we may assume fij = 1. All the 
file's and fij's which have not been fixed at this point are then chosen randomly 
and independently. All coefficients are to be elements in F q . If a solution to 
the network connection problem exists with the a^j's and the fij's specified, 
it is possible to determine values of &t,j,j at the sinks such that a solution to 
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the network connection problem is given. Let fj, be the number of variables 
djj and fij chosen randomly. Call these variables X\, . . . ,X^. Consider the 
polynomial P in ([2]) and let P be the polynomial made from P by plugging in 
the fixed values of the a^/s and the fixed values of the fij's (calculations taking 
place in F g ). Then P is a polynomial in X±, . . . , X^. The coefficients of P are 
polynomials in the 6t.jj's over F q . Finally, define 

P :=P rem (Xf -X U ...,X*- XJ. 

The success probability of random linear network coding is the probability that 
the random choice of coefficients will lead to a solution of the network connection 
problerrQ as in Section [2] That is, the probability is the number 



|{(x 1 ,...,^)GF^|P(x 1 ,...,^)^0}| 

HI 

(3) 



\{(x 1 ,...,x fl )eF q i \P(x 1 ,...,x IJ ,)=/=0}\ 



qi> 

To see the first result observe that for fixed (xi, . . . , x M ) £ F^ 1 , P(xi, . . . , x M ) 
can be viewed as a polynomial in the variables bt^j's with coefficients in F q and 
recall that the &t,i,j's occur in powers of at most 1. Therefore, if P(xi, . . . , x M ) ^ 
0, then by Proposition [2] it is possible to choose the b^ij's such that if we plug 
them into P(x±, . . . , x M ) then we get nonzero. The last result follows from the 
fact that P(xi, . . . , x M ) — P{x\, . . . , x M ) for all (xi, . . . , x M ) s F^. In this section 
we shall present a method to estimate the success probability using Grobner 
basis theoretical methods. 

We briefly review some basic definitions and results of Grobner bases. See 
[3] for a more detailed exposition. Let Ai(Xi, . . . , X n ) be the set of monomials 
in the variables X\ , . . . , X n . A monomial ordering -< is a total ordering on 
M(Xi, . . . , X n ) such that 

L <M => LN -< MN 



holds for all monomials L, M, N 6 A4(Xi, . . . , X n ) and such that every nonempty 
subset of A4(Xi, . . . , X n ) has a unique smallest element with respect to -<. The 
leading monomial of a polynomial F with respect to -<, denoted by lm(F), is 
the largest monomial in the support of F. Given a polynomial ideal I and a 
monomial ordering the footprint A^(J) is the set of monomials that cannot 
be found as leading monomials of any polynomial in /. The following proposi- 
tion explains our interest in the footprint (for a proof of the proposition see [21 
Pro. 8.32]). 

Proposition 3. Let F be a field and consider the polynomials Fi,...,F g € 
F[Xi, . . . , X n ]. Let L = (Fi, . . . , F s ) C F[Xi, . . . , X n ] be the ideal generated by 
Fi, ... ,F S . If A_<(/) is finite then the number of common zeros of Fi, . . . , F s in 
the algebraic closure ofF is at most equal to |A^(/)|. 

Proposition's known as the footprint bound. It has the following corollary. 

lr This corresponds to saying that each sink can recover the data at the maximum rate 
promised by network coding. 
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Corollary 4. Let F G F[Xi, . . . , X n ] where F is a field containing F q . Fix a 
monomial ordering and let 

Xf ■■■Xt = lm(F rem (Xf -X U ...,X« - X n )). 

The number of zeros of F over F q is at most equal to 

n 

g B - 11(9 (4) 

Proof. We have 

A^((F,X?-X 1 ,...,X«-X n )) 

C A^((lm(F rem (Xf - X u . . . , X* - X n )), Xf, . . . , X«» 

and the size of the latter set equals (|4]). The result now follows immediately 
from Proposition [3] □ 

Theorem 5. Let as above P be found by plugging into P some fixed values 
for the variables a%j, (i,j) € K, and by plugging into P some fixed values for 
the variables fij, G J, and by leaving the remaining /! variables flexible. 

Assume as above that there exists a solution to the network connection problem 
with the same values for these fixed coefficients. Denote by X\, . . . , X^ the 
variables to be chosen by random and define P := P rem (Xf— X%, . . . , X^—X^). 
(Note that if q > \T\ then P — P). Consider P as a polynomial in the variables 
Xi, . . . , X^ and let ~< be any fixed monomial ordering. Writing Xf 1 ■ ■ ■ Xjf = 
LM(P) the success probability is at least 

q-»H{q-3v)- (5) 

v=l 

As a consequence the success probability is in particular at least 



q p min < ]^[(<Z — Sj) Xf 1 ■ ■ • X^ is a monomial in the support of P>. 
^ i=i ' 



(6) 



Proof. Let F be the quotient field F g (Xi, . . . , X^). The result in (O now follows 
by applying Corollary 2] and ([3]) . As the leading monomial of P is of course a 
monomial in the support of P ([5]) is smaller or equal to (JSJ. □ 

Remark 6. The condition in Theorem [5| that there exists a solution to the 
network connection problem with the coefficients corresponding to K and J being 
as specified is equivalent to the condition that P =/= 0. 

We conclude this section by mentioning without a proof that Grobner basis 
theory tells us that the true success probability can be calculated as 

- |A^((P, Xf - Xi, . . . , X« - X,})\). 

This observation is however of little value as it seems very difficult to compute 
the footprint 

A H ({P,Xf-X 1 ,...,X«-X (i }) 
due to the fact that \i is typically a very high number. 
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5 The Bound by Ho et al. 



In [7J Ho et al. gave a lower bound on the success probability in terms of the 
number of edges j with associated random coefficient^ {a^, fi,j}- Letting n be 
the number of such edges [7J Th. 2] tells us that if q > \T\ and if there exists a 
solution to the network connection problem with the same values for the fixed 
coefficients, then the success probability is at least 

(7, 



The proof in [7J of ([7]) relies on two lemmas of which we only state the first one. 

Lemma 7. Let n be defined as above. The determinant polynomial of M t has 
maximum degree r\ in the random variables {cLi,j,fi,j} an d * s linear in each of 
these variables. 

Proof. See [7J Lem. 3] . Alternatively the proof can be derived as a consequence 
of Theorem fTTI in Section 171 □ 

Recall, that the polynomial P in ([2]) is the product of the determinants 
|M t |, t G T. Lemma [7j therefore implies that the polynomial P has at most 
total degree equal to |T|?7 and that no variable appears in powers of more than 
\T\. The assumption q > \T\ implies P — P which makes it particular easy to 
see that the same of course holds for P. Combining this observation with the 
following lemma shows that the numbers in © and © are both at least as large 
as the number [[7]). 

Lemma 8. Let 77, \T\,q <E N, |T| < q be some fixed numbers. Let fx, x±, . . . , x^ G 
No satisfy 

0<xx<\T\,...,0<x^<\T\ 
and Xi + ■ ■ ■ + Xfj, < | T 1 77 . The minimal value of 

n / 

q-Xi 



n 



(taken over all possible values of n, x±, . . . , x^) is 

\T\ x 



Proof. Assume n and xy , . . . , x^ are chosen such that the expression attains its 
minimal value. Without loss of generality we may assume that 

xi > x 2 > ■ ■ ■ > x^ 

holds. Clearly, x\ + • • • + x^ = \T\n must hold. If Xi < \T\ and Xj+i > then 

(q - Xi)(q - x i+1 ) > (q - (Xi + l))(q - (x i+1 - 1)) 

which cannot be the case. So x\ = ■ ■ • = x^ = \T\. The remaining Xj's if any 
all equal zero. □ 



2 We state Ho et al.'s bound only in the case of delay-free acyclic networks. 
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Figure 1: The network from Example [5] 

6 Examples 

In this section we apply the methods from the previous sections to two concrete 
networks. We will see that the estimate on the success probability of random 
linear network coding that was described in Theorem O can be considerably 
better than the estimate described in [7J Th. 2] . Also we will apply the method 
from Section [3] to determine the smallest field of characteristic two for which 
network coding can be successful. 

As random linear network coding is assumed to take place at the nodes in a 
decentralized manner, one natural choice is to set fij = 1 whenever the indegree 
of the end node of edge i is one and j is the downstream edge adjacent to i. 
Clearly, if j is not a downstream edge adjacent to i we set /jj = 0. Whenever 
none of the above is the case we may choose fij randomly. Also if there is only 
one source and the outdegree of the source is equal to the number of symbols to 
be send we may enumerate the edges from the source by the numbers 1, . . . , h 
and set ajj = lifl<i = j<h and set = otherwise. This strategy 
can be generalized also to deal with the case of more sources. In the following 
two examples we will choose the variables in the manner just described. The 
network in the first example is taken from [H Ex. 3.1] whereas the network in 
the second example is new. 

Example 9. Consider the delay-free and acyclic network in Figure[]i There 
is one sender v\ and two receivers v\2 and V13. The min-cut max- flow number 
is two for both receivers so we assume that two independent random processes 
emerge from sender v\. We consider in this example only fields of characteristic 
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Table 1: From Example [5J Estimates on the success probability 



q 


4 


8 


16 


32 


64 


-Pnew l(<?) 


0.140 


0.430 


0.672 


0.893 


0.909 


Pnew 2 (q) 


0.703 x 10- 1 


0.322 


0.588 


0.773 


0.880 


Pno(q) 


0.625 x 10" 1 


0.316 


0.586 


0.772 


0.880 



2. Following the description preceding the example we set a^i = 0,2.2 = 1 and 
aij = in all other cases. Also we let fij = 1 except 

f3,7, h,7, /4,10, /g.lOi J9,ll! J6,ll! /l2,16j /l4,16 

which we choose by random. As in the previous sections we consider bt.i,j a>s 
fixed but unknown to us. The determinant polynomial becomes 

2„2„2„i, 1 „2j2„l 1 „2j2j-2. 



P = {b 2 c 2 e 2 gh + c 2 f 2 gh + a 2 d 2 f 2 gh)Q 



where 



a = /3,7 & = /5,7 c = / 4 ,io d = /s,10 

e = /g,n / = /6,n .9 = Zi2,i6 h = fu,ie 

and Q = \B V \ \B' Vl3 \. Here, B' Vl2 respectively B' v is the matrix consisting 
of the nonzero columns of B Vl2 respectively the nonzero columns of B Vli . Re- 
stricting to fields F q of size at least 4 we have P — P and we can therefore 
immediately apply the bounds in Theorem Applying © we get the following 
lower bound on the success probability 

(q-2f{q-lf 
inew 2{q — c • 

g 5 

Choosing as monomial ordering the lexicographic ordering -<i ex with 
a -<iex b ^icx d -Kiex e -<i cx 9 -<iex h -<i cx f -<i cx c 

the leading monomial of P becomes c 2 f 2 gh and therefore from ([5]) we get the 
following lower bound on the success probability 

(q-2f{q-lf 
r new 1\M} — a 

q 4 

For comparison the bound ([7|) from jTjj states that the success probability is at 
least 

We see that P new 1 exceeds Pho with a factor (q — l) 2 /{q — 2) 2 , which is larger 
than 1. Also P new 2 exceeds Pro- In Table [7] we list values of P new i(q), 
Pnew s{q) and P# (g) for various choices of q. 

We next consider the field F2. We reduce P modulo (a 2 — a, . . . ,h 2 ~ h) to 

get 

P = (bcegh + cfgh + adfgh)Q. 
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Figure 2: The network from Example [TU] 



From (jnj) we see that the success probability of random network coding is at 
least 2 -5 . Choosing as monomial ordering the lexicographic ordering described 
above ([5j) tells us that the success probability is at least 2~ 4 . For comparison 
the bound ([7]) does not apply as we do not have q > \T\. It should be mentioned 
that for delay-free acyclic networks the network coding problem is solvable for 
all choices of q > \T\ J3J/ and Jllf . From this fact one can only conclude that 
the success probability is at least 2~ 8 (8 being the number of coefficients to be 
chosen by random). 

Example 10. Consider the network in Figured The sender v\ generates 3 
independent random processes. The vertices V\\, Vn and V\i are the receivers. 
We will apply network coding over various fields of characteristic two. We start 
by considering random linear network coding over fields of size at least 4. As 
4 > |T| =3 we know that this can be done successfully. 

We set ai^i = 02,2 = &3.3 = 1 and a^.j = in all other cases. We let = 1 
except / 4 ,i 3 , / 7j i3, / S) i4, / 8) i4, /i ,i4, 79,15, /n,i5, which we choose by random. As 
in the last section we consider &t » j as fixed but unknown to us. Therefore P = P 
is a polynomial in the seven variables / 4 ,i3, /7,13, /s,i4> /s.14, /io,U, /9,1s, /n,i5- 
The determinant polynomial becomes 

P = (abcdefg + abce 2 f + b 2 c 2 efg)Q, 

where 

a = /4,13 b = / 5j i4 c = / 7i i3 d = f 8M 

e = /g,i5 / = /1044 g — /1145 

and Q = \B' Vii \ \B' Vi2 \ \B' Vl3 \. Here, B' V11 respectively B' Vi2 respectively B' Vi3 is 
the matrix consisting of the nonzero columns of B Vll respectively the nonzero 
columns of B Vl2 respectively the nonzero columns of B Vl3 . Choosing a lexico- 
graphic ordering with d being larger than the other variables and applying ([5]) 
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Table 2: From Example [TU1 Estimates on the success probability 
q 4 8 16 32 64 

P now i(q) 0.133 0.392 0.636 0.800 0.895 
fncw 2 (<7) 0.105 0.376 0.630 0.799 0.895 
PBo(q) 0.156 x 10- 1 0.244 0.536 0.744 0.865 



we get that the success probability is at least 

,_ (g-l) 7 

*new 1\Q) 

Applying ([HJ) we see that the success probability is at least 

(q - l) 3 (g - 2) 2 



Pnew 2(1) — 



For comparison tells us that success probability is at least 

(9-3) 3 



P Ho {q) 



<? 3 



Both bound §5§ and bound ^ exceed ([7]) for all values of q > 4. In Table\Mwe 
list Pnew i(q)> Pnew 2(9) and Ph (<i) for various values of q. 

We next consider the field F2. We reduce P modulo (a 2 — a, . . . , g 2 — g) to 

get 

P = (abedefg + abcef + bcefg)Q. 

From ^ we see that the success probability of random network coding is at least 
2 -7 . Choosing a proper monomial ordering we get from ^ that the success 
probability is at least 2 -5 . For comparison neither J^, nor Ml] tells us that 
linear network coding is possible. 

7 The Topological Meaning of \M t \ 

Recall from Section [5] that Ho et al.'s bound ([7]) relies on the rather rough 
Lemma [7j The following theorem gives a much more precise description of 
which monomials can occur in the support of P and P by explaining exactly 
which monomials can occur in \M t \. Thereby the theorem gives some insight 
into when the bounds ([5|) and are much better than the bound (|7|). The 
theorem states that if if is a monomial in the support of |M t | then it is the 
product of flip's, /ij's and &t,ij's related to h edge disjoint paths P\,...,Ph 
that originate in the senders and end in receiver t. 

Theorem 11. Consider a delay- free acyclic network. If K is a monomial in 
the support of the determinant of M t then it is of the form K\ ■ ■ ■ Kh where 

R v = //(">, 4") //<"),;<") ' ' ' ^il > _ 1 ,iil > Kv u ,i[l > 

for u — 1, . . . , h. Here, {v±, . . . , v^} = {1, . . . , h} holds and l± , . . . , zi respec- 
tively zli , . . . , li^ are pairwise different. Further 
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Figure 3: The butterfly network 



unless u\ — ui and i — j hold. In other words K corresponds to a product of h 
edge disjoint paths. 



Proof. A proof can be found in the appendix. 
We illustrate the theorem with an example. 



□ 



Example 12. Consider the butterfly network in Figure [3J A monomial K is 
in the support of \M tl | if and only if it is in the support of the determinant of 



N t! = (ni,j) 



I + F Bl 
A 



1 





/l,3 


flA 


























1 








/2,5 


/2,6 























1 




















& ii,l,3 


6*1,2,3 











1 








/4,7 


























1 





/5,7 





























1 



































1 


/7,8 


JV,9 





























1 





6*1,1,8 


6*1.2,8 


























1 








01,1 


Ol,2 





























02,1 


02,2 






























By inspection we see that the monomial 

K = 01,102, 2^*i, 1,3^*1, 2,8^7, 8.f5, 7/2, 5.fl, 3 

is in the support of \N tl \ . We can write K = K\K% where 

^1=01,1/1,36(^1,3 and A" 2 = 02,2/2, 5/5, 7/7,86*1, 2,8- 

This is the description guaranteed by Theorem To make it easier for the 
reader to follow the proof of Theorem \ll\ in the appendix we now introduce some 
of the notations to be used there. By inspection the monomial K can be written 



K 



11 

11 n 

i=l 



i,p{i) 



12 



where the permutation p is given by 



p(l) = 3 p(2) = 5 p(3) = 10 p(4) = 4 p(5) = 7 p(6) = 6 
P(7) = 8 p(8) = ll p(9)=9 p(10) = l p(ll)=2 

Therefore if we index the elements in {1, . . . , 11} &?/ 

ii = 10 %2 = 1 «3 = 3 «4 = 11 «5 = 2 «6 = 5 

i? = 7 is — 8 ig = 4 iio = 6 «n = 9 
t/ien we can write 

K\ n ii ,p(ii) n i2 ,p{i2) n i3 ,p(i3) 

K2 — n i 4 , p(i 4 ) n i 5 , p(i 5 ) n i e , p(i 6 ) n i 7 ,p(i 7 ) n is,p(i 8 ) 

and we have 

n i9,p(ia) ~ n iia,p{iia) ~ n in,p(in) — ^ 

corresponding to the fact p(ig) = £g, p(iio) = *io and p(in) = in. 

Remark 13. T/ie procedures described in the proof of Theorem \11\ can be re- 
versed. This implies that there is a bijective map between the set of edge disjoint 
paths Pi, . . . ,Ph in Theorem ] 11\ and the set of monomials in \M t \. 

Theorem 1111 immediately applies to the situation of random network coding 
if we plug into the a,j's and into the /t.ij 's on the paths Pi, . . . ,Ph the fixed 
values wherever such are given. Let as in Lemma [7] 77 be the number of edges 
for which some coefficients djj , fi j are to be chosen by random. Considering 
the determinant as a polynomial in the variables to be chosen by random with 
coefficients in the field of rational expressions in the &t,j,j's we see that no mono- 
mial can contain more than 77 variables and that no variable occurs more than 
once. This is because the paths Pi, . . . , Ph are edge disjoint. Hence, Lemma [7] 
is a consequence of Theorem QTJ 

8 Acknowledgments 

The authors would like to thank the anonymous referees for their helpful sug- 
gestions. 



A Proof of Theorem [TT1 

The proof of Theorem [Tl] calls for the following technical lemma. 

Lemma 14. Consider a delay-free acyclic network with corresponding matrix 
F as in Section^ Let I be the \E\ x \E\ identity matrix and define 

r = ( 7iJ ) = j + F. 

Given a permutation p on {1, . . . ,\E\} write 

i times 

P W (A)=p(p(---(A)---)) 
If for some A £ {1, . . • , \E\} the following hold 
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(1) X,p{\), . . . ,p( x '(\) are pairwise different 

(2) p(-+ 1 )(A)e{A,p(A),...,p(-)(A)} 



(3) 7A,p(A),7p(A) : p(p(A)), 

then x = 0. 



>7 P (-)(A),p<-+i)(A) are a ^ nonzero 



Proof. Let p be a permutation and let x and A be numbers such that (1) (2) 



and |(3)| hold. As p is a permutation then (1) and (2) implies that (A)) = A 

Aiming for a contradiction assume x > 0. As p(n) — r\ does not hold for any 

v e{\,p(\),..., P (*\\)}, 



7A,p(A) i 7p(A),p( 2 ) (A) i ■ ■ ■ i 7p(*) (A),p(^ + 1 > (A) 

are all non-diagonal elements in I + F. By |(3)| we therefore have constructed a 
cycle in a cycle-free graph and the assumption x > cannot be true. □ 

of Theorem \lll A monomial is in the support of the determinant of M t if and 
only if it is in the support of the determinant of 

A7 (I + F Bj\ , , 

To ease the notation in the present proof we consider the latter matrix. Let p 
be a permutation on {1, ... , \E\ + h} such that 

\E\+h 

n n ^(s) + °- ( 8 ) 

s=l 

Below we order the elements in {1, . . . , \E\ + ti} in a particular way by indexing 
them ix, . . . , i\E\+h according to the following set of procedures. 
Let ix = \E\ + 1 and define recursively 

i s = p(i s -i) 

until \E\ < p(i s ) < \E\ + h. Note that this must eventually happen due to 
Lemma IT4"1 Let si be the (smallest) number such that \E\ < p(i Sl ) < \E\ + h 
holds. This corresponds to saying that n^^fa) is an entry in A, that «i 2 . p (i 2 ), 
■ ■ ■ j ri i si _i.p(i si _i) a re entries in / + P, and that ni siiP (i si ) is an entry in Pyf . 
Observe, that p(i r ) = i r cannot happen for 2 < r < s\ as already p(i r -i) = i r 
holds. As n ir lP u r ) is non-zero by ([5]) we therefore must have 

^i r ,p(i r ) fi T ,p(i r ) Ji r ,i r _)-i 

for 2 < r < sx- Hence, 

(^ii ,p(ii) i ■ • ■ ? n i si ,p(i si ) ) = ( a M2 ' fi2,is > * * • ' fis 1 -l,is 1 > ^t,Vi ,i sl ) 

for some i>x- Denote this sequence by Px. Clearly, Px corresponds to the poly- 
nomial in the theorem. 

We next apply the same procedure as above starting with i Sl +i = \E\ + 2 
to get a sequence Pi of length S2- Then we do the same with ? Sl+S2+ i = \E\ + 
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3, . . . , i Sl+ ... Sh _ 1+ i = \E\ + h to get the sequences P3, . . . , Ph- For u = 2, . . . , h 
we have 

^ u = ( n ia 1 + .-. + Su _ 1 + l,p(i sl +- - + Su _ 1 +l)' • • • ! n i, 1 + ... +Su ,p(is 1+ .-. + s u )J 

= ^tl,t, 1 + ... + >ii _ 1+ 2, /i sl + ... +au _ 1+2 ,lsi + --- + s tt _ 1 + 3! ■ • • 1 

/f«i +— +*tt-ij**i+---+»« ' bt,v u ,i sl + ... + . 

Clearly, P u corresponds to K u in the theorem. Note that the sequences Pi , . . . , Ph 
by the very definition of a permutation are edge disjoint in the sense that 

(1) riij occurs at most once in Pi, . . . , Ph, 

(2) if rij^, rijj 2 occur in Pi, . . . , Ph then l\ = I2, 

(3) if n jl . h n h .i occur in Pi, . . . ,P h then ji = j 2 . 

Having indexed si + • • • + Sh of the integers in {1, . . . , \E\ + h} we consider what 
is left, namely 

A = {l,...,\E\ + h}\{i 1 ,...,i Sl+ „. +Sh }. 

By construction we have i\ — \E\ + 1, . . . , i Sl _| |-s h _i+i — |P| + h and therefore 

AC{1,...,|£?|}. Also by construction for every 

S e {l,...,\E\}n{h,...,i Sl+ ... +8h } 

we have 5 = p(e) for some e € . . . , i si _| ^ Sh }- Therefore p(A) £ A for all 

A G A holds. In particular p^ x '(X) € {1, . . . , \E\} for all x. From Lemma IT4l we 
conclude that p(X) = A for all A 6 A. □ 
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